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Nuclear modification factor (Rep) of protons and pions are investigated by simulating An + Au 
collisions from 0.8 to 1.8A GeV in a framework of an isospin-dependent quantum molecular dynamics 
(IQMD) model. Rep of protons rises with the increase of pr at different beam energies owing to 
radial flow and Cronin effect. The rate of increase of Rep is suppressed at higher beam energies. The 
significant difference of Rep between protons and pions indicates different medium effects between 
protons and pions. By changing the in-medium nucleon-nucleon cross section, the Rep of protons 
changes a lot, while the Rep of pions does not. Taking the pion absorption into account, the Rep of 
pions becomes close to unity without pr dependence after deactivating the reaction nN —>■ A, while 
there is nearly no change on proton. This suggests that the pion absorption plays a dominant role 
on pion dynamics and have slight effect for proton dynamics. 

PACS numbers: 25.70.-z, 21.65.Mn 


I. INTRODUCTION 

One of the main goals of research in the intermediate 
energy heavy ion collisions (HIGs) has focused on learn¬ 
ing the bulk properties of hot and compressed nuclear 
matter and its transport mechanisms since the last thirty 
years [ll-[l|- Transport models such as BUU-type [l-@ 
and Molecular Dynamics type 0,101 have been suc¬ 
cessful in describing the reaction dynamics of the low and 
intermediate energy heavy-ion collisions. The two main 
ingredients of the nuclear transport process are the nucle¬ 
onic mean field and nucleon-nucleon (NN) binary inter¬ 
action. Recently, the medium effects on nucleon-nucleon 
cross section (NNGS) have been widely investigated by 
replacing the NNGS in vacuum with an in-medium one 
Since a high density region for the compressed 
nuclear matter could be reached up to 2-3 times nor¬ 
mal nuclear matter density po before it expands during 
the process of heavy-ion collision at 1-2A GeV, the in¬ 
medium NNGS is therefore an important component in 
these phenomenological simulations due to its close rela¬ 
tion with the density. 

The nuclear modification factor (NMF) has been ex¬ 
tensively studied in relativistic heavy ion collisions in re¬ 
cent years [Hill. In these studies, unanimous results 
have demonstrated that the NMF is suppressed at high 
Pt owing to partonic energy loss effect. However, the Rep 
of protons show a rise with pp at the low and moderate pp 
range in intermediate energy HIGs, which was argued to 
be an indication of combined effect from radial flow and 
Gronin effect [l^. Nuclear modification factor of light 
nuclei has been also investigated for the first time in a 
framework of thermal and coalescence models and a num- 
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ber of constituent quark scaling like behavior of pions and 
protons was exhibited [2^. On the other hand, the nu¬ 
clear stopping can provide the information on the nuclear 
equation of state (EoS), in medium nucleon-nucleon cross 
section as well as the degree of equilibrium [2^ . Fur¬ 

thermore, the magnitude of nuclear stopping may have 
a direct relationship with the enhancement of Rep- Be¬ 
sides, some other physical quantities, such as radial flow, 
temperature and viscosity, can also provide abundant in¬ 
formation about dense hadronic matter formed in heavy 
ion collisions [26l - l2^ . 

In the present work, the nuclear modification factors 
Rep of protons and pions at different incident energies are 
investigated systematically. The nuclear medium effect 
from in-medium NNGS has been studied, while the pion 
absorption effect has also been discussed. 

The article is organized as following: In Sec.II a brief 
introduction on IQMD model is described. Sec. Ill gives 
a description of pion dynamics in the IQMD model from 
which the pion spectra are compared with the FOPI 
experimental results and a preferable matching is ob¬ 
tained. The Rep of protons and pions versus pp and 
the radial flow are calculated at different incident energy 
in Sec. IV. Then we turn to the study of in-medium 
nucleon-nucleon cross section which indicates the effec¬ 
tive nucleon-nucleon interaction in the dense matter en¬ 
vironment in Sec. V. At last, the effect of pion absorption 
has been investigated in Sec. VI. Summary is given in 
the last section. 


II. BRIEF DESCRIPTION OF IQMD MODEL 

The Quantum Molecular Dynamics model is a trans¬ 
port model which is based on a many body theory to 
describe heavy ion collisions from intermediate to rela¬ 
tivistic energy 01113, mi- An extended version, so- 
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called the Isospin dependent Quantum Molecular model 
(IQMD) which considers the isospin effects is suitable 
to investigate asymmetric nuclear system. IQMD can 
successfully treat collective flow, multi-fragmentation, 
isospin effects, transport coefficients, giant resonance, 
and strangeness production etc [H, . 

Wave function of each nucleon in the IQMD model is 
described as a coherent state with the form of Gaussian 
wave packet. 


Mr.t) = 


(27rL)3/4 


expi -- )exp{ - - -), 


a) 

where ri and pi are the time dependent variables which 
describe the center of the packet in coordinate and mo¬ 
mentum space, respectively. The parameter L, related to 
the width of wave packet in coordinate space, is deter¬ 
mined by the size of reaction system. Usually L = 2.16 
fm^ for Au-|-Au system. The wave function of the sys¬ 
tem is the direct product of all the nucleon wave functions 
without considering the Fermion property of nucleon: 


action potential in Eq. |4]can be expressed as follows: 
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The symmetry potential between protons and neutrons 
corresponding to the Bethe-Weizsacker mass formula can 
be taken as 


^sym ~ ^^T'4iT^j5{ri Tj) (7) 

with t6=100 MeV. By integrating Skyrme part as well as 
the momentum dependent part of the two-body interac¬ 
tion and introducing the interaction density. 


^r,t) = l[Ur,t). ( 2 ) 
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As a compensation, Pauli blocking is employed in the 
initializations and collision process to restore some parts 
of the quantum property of many Fermion system. 

By applying a generalized variational principle on the 
action of the many-body system, one can get the equa¬ 
tions of motion for pi and ri, which are listed as follows 




dn 


d{H) 

dpi 


( 3 ) 


The Hamiltonian {H) = (T) -|- (U) where T is the ki¬ 
netic energy, the potential V is expressed by 


one can get the local mean field potential which contains 
the Skyrme potential and momentum dependent poten¬ 
tial 


U = — + — [ dpg{p')5\n[e{p-p')^ + l], 

Po Po Po J 

(9) 

where po is the saturation density at ground state, 
g{p,t) = e ’ is the momentum dis¬ 

tribution function, the interaction density p = '^ij Pij , 
and a, /3, and 7 are the Skyrme parameters, which con¬ 
nect tightly with the EOS of the bulk nuclear matter, as 
listed in Table HI 


Mr,p,t) fj{r',p',t) df dr’dpdp', 

i j^i 

( 4 ) 

where the Wigner distribution function fi(f,p,t), which 
is the phase-space density of the ith nucleon, is obtained 
by applying the Wigner transformation on the single nu¬ 
cleon wave function: 


fi{f,p,t) = 


{irhy 




( 5 ) 


The baryon-potential consists of the real part of the 
G-Matrix which is supplemented by the Coulomb inter¬ 
action between the charged particles. The former one can 
be divided into three parts, the Skyrme-type interaction, 
the finite-range Yukawa potential, and the momentum- 
dependent interaction (MDI) parts. The two-body inter- 


TABLE I: Parameter sets for the nuclear equation of state 
used in the IQMD model. S and H represent the soft and 
hard equation of state, respectively, M refers to the inclusion 
of momentum dependent interaction. This table is adapted 
from [1]. 



a 

P 

7 


e 


(MeV) 

(MeV) 


(MeV) 

( (GeV\-‘ ) 

s 

-356 

303 

1.17 

- 

- 

SM 

-319 

320 

1.14 

1.57 

500 

H 

-124 

71 

2.00 

- 

- 

HM 

-130 

59 

2.09 

1.57 

500 


With the help of coalescence mechanism, the informa¬ 
tion of fragments produced in HICs can be identified in 
IQMD. A simple coalescence rule to form a fragment is 
used with the criteria Ar = 3.5 fm and Ap = 300 MeV/c 
between two considered nucleons. 
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III. DELTA RESONANCE PRODUCTION AND 
DECAY, AND PION DYNAMICS 

In the present work, the production of pion is consid¬ 
ered when inelastic scattering occurs. Pions are produced 
via the decay of A resonance, the following inelastic reac¬ 
tion channels have been taken into account explicitly at 
lA GeV domain [s^. And the higher mass resonance N* 
channel hasn’t been considered owing to its contribution 
is negligible in this energy domain. 

{a)NN —>■ NA (hard A production), 

(6 ) A -)■ Att (A decay), , . 

(c)AN —>■ NN (A absorpion), ^ ' 

{d)NTr —>■ A (soft A production). 

Elastic TT — TT, tt—N, tt—A, A —A and A — N scattering 
channels are not taken into account. In the processes (a) 
and (d), the experimental cross section and the elastic 
NN collision are used as shown below. 

The elastic nucleon-nucleon scattering angular distri¬ 
bution is described as ^ e'^ *, where t is equal to 
—2p^(l — cos(0)), p is the total momentum in the center 
of mass system (c.m.s) of these two colliding nucleons, 
the squared momentum transfer c is a function of ^/s 
and defined as, 


(3.65-(^- 0.18766))^ 
UVA c (l + (3.65-(yi- 0.18766))6)’ 


( 11 ) 


where is the c.m.s energy with GeV unit [3^ . 

The inelastic angular distribution is given by the fol¬ 
lowing function 


dCel 

dfl 




( 12 ) 


where the parameters a(^/s) and b(y/s) were suggested 
by Huber and Aichelin and vary in their definition for 
different intervals of ^/s. The details can be found in 
Ref. [H]. 

Pions produced via A decay propagate with a high 
thermal velocity under the Coulomb force. The different 
isospin channels have been considered and the branching 
ratios use the Clebsch-Gorden coefficient: 

(a) A++—>■ l(p-I-7r+) 

(b) A+—> 2/3(p-I-7r°)-I-l/3(n-I-7r+) 

(c) A° —2/3(n -I-7 r°) -I- l/3(p -I-7 r“) (13) 

(d) A^ —>■ l(n-I-7r“) 


the TT spectra from the IQMD model calculation by a 
soft equation of state with momentum dependent inter¬ 
action (SM-EOS) are compared with the EOPI result: 
it shows the yields of charged and neutral tt measured 
around mid-rapidity (±0.2) in Au-b Au collisions at lA 
GeV by MMD simulation are comparable with the EOPI 
results |39| . Under the “minimum bias” condition, pi¬ 
ons are chosen in the rapidity range —0.2 < Y < 0.2, 
and normalize to dy=l. Figure [T] shows the spectra can 
match well with experimental data except a little under¬ 
estimated in higher pT region because all pions are only 
produced by A decay in our QMD model calculations. 



FIG. 1: (Color online) Comparison of the invariant tt pro¬ 
duction cross section of the measurement and our simulations 
for mid-rapidity pions in Au + Au collisions at lA GeV. The 
circle represents the FOPI data of tt~ and the square for 
the red, green and blue lines represent different types of tt 
from our IQMD simulations. 


IV. NUCLEAR MODIFICATION FACTOR R^p 
OF PROTONS AND PIONS 


To study nuclear medium effect in Au-|-Au collisions, 
it is convenient to introduce a ratio Rep of the particle 
yiel d in central collisions to that in peripheral collisions 
(la liol . Both yields are normalized by corresponding 
nucleon-nucleon binary collision numbers (Ncoii) (binary 
scaling), 


After a TT is produced, it has two routines, namely tt 
absorption ttNN — >• AN — >■ NN and tt scattering (reab¬ 
sorption) ttN —>■ a —>■ ttN. In this calculation, we named 
the sum of this two routines “tt absorption process”. 

In our previous publication, the px spectra of light 
charged particles (p, d, t) have been shown to repro¬ 
duce the experimental data very well (2^. In this work. 


Yield(central)/Nf^n 
Yield(peripheral) / N^^n 

If nucleus-nucleus collision is a mere superposition of 
Neou independent nucleon-nucleon collisions, the R^p 
would be equal to unity. However, owing to the medium 
effect, such as initial multiple scattering, secondary par- 
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tide decay and shadowing effect etc, the NMF will show 
deviation from unity. 

In the present work, Au+Au collisions at 0.8, 1.0, 1.2, 
1.5 and 1.8A GeV are simulated with IQMD model for 
the soft equation state with momentum dependent in¬ 
teraction (SM-I-MDI). The double differential transverse 
momentum spectra [ 2 -w^dpTdy ) protons (pions) at dif¬ 
ferent centralities (0 — 20%, 20 — 40%, 40 — 80%) have been 
obtained with the c.m.s. rapidity cut {\Y/Yproj\ < 0.1, 
where Yp^oj is the initial projectile rapidity). All the Rep 
in this work are obtained by dividing the spectra in the 
centrality of 0-20% to the one in the centrality of 40-80%. 


A. Rep of protons and pions at different energies 

Figure [2] shows the Rep of protons and pions for en¬ 
ergies 0.8-1.8 AGeV. It is seen that the Rep of protons 
enhances quickly as px increases at different beam ener¬ 
gies in Fig. [2Ia), while the Rep of tt increases at low px 
and level off at high px as shown in Figl^Jb). For pro¬ 
tons, the strength of Rep enhancement is suppressed at 
high Px with the increasing of beam energy. We noticed 
that this energy dependence of Rep is quite similar to 
the preliminary results of Rep [4l| which show a mono¬ 
tonic evolution with collision energy from enhancement 
at low energy due to the Gronin effect to suppression at 
high energy due to the partonic energy loss when jets 
pass through the hot dense quark-gluon matter. On one 
hand, in the low px region, radial flow plays a major 
role in central collisions, which pushes protons to higher 
Px region and results in the smaller Rep at low px- On 
the other hand, in the high px region, the Gronin effect 
due to multiple nucleon-nucleon scattering effect 
tends to transform the longitudinal momentum into the 
transverse momentum and the effect becomes stronger 
with the increasing of px in central HIGs, which leads 
to larger Rep at high px- For the pion case, however, 
the Rep of pions has very different dynamical origin be¬ 
cause they are produced via hard A decay in NN collision 
and absorbed by nucleon with a large probability (81% 
in Au-|-Au at lAGeV), [l^. These differences mean the 
different nuclear medium effects for protons and pions. 


B. Radial flow 

Two important physical quantities, nuclear stopping 
and the radial flow, which are sensitive to the proper¬ 
ties of nuclear bulk matter, have been extensively inves¬ 
tigated in a wide range of incident energy from tens of 
MeV/nucleon to hundreds of GeV/nucleon in many ex¬ 
perimental or theoretical work [2^, IdSl - ld^ . The nuclear 
stopping can be described with a ratio of transverse to 
parallel quantities (energy or momentum), it reflects how 
much energy of original longitudinal motion is transferred 
into the internal degrees of freedom. 




FIG. 2: (Color online) Rep of protons (a) and pions (b) versus 
Pt at different beam energies. Shadow in lines represent the 
statistical error. 


In the compression stage, with nucleon-nucleon r fre¬ 
quently occurring, a highly dense and thermal nuclear 
matter zone is formed. New species of particles (such as 
tt) can be created in the bulk matter with unique condi¬ 
tion, and their production yield or emitting pattern can 
be used to probe the global properties of their surround¬ 
ing, just as the Rep of pion in the previous Sec. HVT Aj. 
and pion absorption effect in Sec. ED 

And then, the expansion stage occurs owing to the 
high pressure in the compressed region. While the parti¬ 
cles on the surface of coupling matter zone emit outward, 
and inner particles frequently interact with each other, 
which cause that the probability of outward motion is 
larger than the one of inward motion, and then, the col¬ 
lective motion of radial flow grows until the nuclear sys¬ 
tem freezes out. At high incident energy, hydrodynamics 
might be suitable to describe these characteristics. A pi¬ 
oneering theoretical model named the Blast-wave model 
has been put forward by Siemens and Rasmussen ^ , 
and lots of work based on it have been carried out 1^, 49 - 


The transverse velocity distribution [3^ in the region of 
0 ^ Rmax is described by a self-similar prohle, which 
is parameterized by the surface velocity /3s: Prip) = 
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FIG. 3: (Color online) Blast-wave fits to the spectra and their 
parameters. Upper pannel: pr spectra in Au-I-Au central 
collisions at 0.8, 1.0, 1.2, 1.5, 1.8A GeV, the solid lines are the 
Blast-wave fits; Lower pannel: the Blast-wave fitting contours 
at different beam energies. 


where Rmax is the freeze-out radius, defined as 
the maximum radius of the expandirng source at ther¬ 
mal freeze-out time, and the Ps is the particles’ radial 
velocity at the maximum surface where radius is equal 
to the freeze-out radius, and the exponent a represents 
the radial flow profile, which describes the evolution of 
the flow velocity with radius (if q;= 0, it means the uni¬ 
form velocity; if a=l, it is similar to Habble’s law; and 
if a=2, it is hydrodynamical evolution). 

Particle spectra are a superposition of individual ther¬ 
mal sources with different r, each boosting with the boost 
angle p = tanh~^Pr{r) [H, 


dn 

prdpT 


oc 



J/ 


rriTCOshp 

Tf 


(15) 

where Ki,Io are the modified Bessel functions and Tf is 
the freeze-out temperature. The shape of spectra is es¬ 
sentially determined by the freeze-out temperature, the 
velocity of the transverse expansion, the flow profile and 
the mass of the particle.The average flow velocity is esti¬ 


mated by taking an average over the transverse geometry: 

iPr) = PS2^- 

Figures [UJa) shows the Blast-wave fitting on the pr 
spectra of the Au-|-Au central collisions at different inci¬ 
dent energies, namely 0.8, 1.0, 1.2, 1.5 and 1.8 A GeV. 
And Fig. 13Kb) are the corresponding fit parameters. From 
the figure, one can see that both the radial flow /3 and the 
freeze-out temperature T, are increase with beam energy. 


V. IN-MEDIUM CROSS SECTION 




FIG. 4: (Color online) Dependence of Rep of protons (a) and 
pions (b) on the reduction factor rj of in-medium NNGS: rj=0 
(black line), r/—0.2 (red line), 77=0.5 (green line) and 77=0.9 
(blue line) 


Usually the free-space nucleon-nucleon cross section 
obtained by experimental measurement is used as 
a default NNCS in QMD model. However, the real in¬ 
medium nucleon-nucleon cross section 
ferent from the free-space nucleon-nucleon cross section 
because of the effects of Pauli blocking and finite system 
of nucleus in heavy-ion reactions etc. In-medium two- 
body cross sections are therefore an indispensable com¬ 
ponent to compensate the nuclear equation of state in the 
QMD simulation [l^, MM- The in-medium nucleon- 
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nucleon cross section can be parameterized from the Par¬ 
ticle Data Group with medium modification which can 
be implemented according to the density dependent pre¬ 
scription [s^: 

<^NN = ■ {1- r]—), (16) 

Po 

where 77 is the in-medium NNCS reduction factor, varied 
between 0 and 1, po is the normal nuclear matter density, 
and p is the local density. Otherwise, the in-medium 
NNCS scaled by the effective mass m*, itnn = ' 

(m *has also been employed in the BUU simulation 
[54| . The latter scaling presumes that, for given relative 
momentum, the matrix elements of interaction are not 
changed between the free space and medium. 



FIG. 6: (Color online) Azimuthal distribution of 7r° with or 
without TT absorption. 
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[^ . In fact, the medium effect is different in various 
ranges of incident energy and matter density (l^ . 


A. Effect on Rep 

Rep is a good quantity to study the effect of the in¬ 
medium NNCS. In Fig. |4l the Rep of protons (upper 
panel) and pions (lower panel) are shown with different p 
values. Due to the limited statistics, there exists fluctua¬ 
tion in high pT region, especially for the 77=0.9 case.The 
Rep of protons has an increasing trend which is explained 
by the Cronin effect . And its trend is increasing more 
rapidly with px in the low 77 value case because of high 
collision rate between nucleons. Collisions become cer¬ 
tainly less in higher 77 values and it makes the Cronin 
effect become less important. On the other hand, the 
trend of pion’s Rep doesn’t seem to have any obvious 
change with different 77 value. The reason might be the 
similar change of cross section for pion interaction in cen¬ 
tral and peripheral collisions when the 77 value changes 
and then lead to the unchanged Rep with 77. However, 
the significant decrease of cross section in the 77 = 0.9 case 
induces a large fluctuation in high px region as shown in 
Fig. mb). In the next section, pion absorption will be 
discussed in detail. 


FIG. 5: (Color online) The Blast-wave fitting to px spectra of 
protons (upper panel) and their fitting parameters of P and 
T (lower panel) with different in-medium NNCS reduction 
factor in the centrality of 0-20%. 

In this part, the aim is to draw a conclusion on the in¬ 
medium NNCS by investigating Rep with the changing 
77 values. The 77 value at 0.2, 0.5, and 0.9 were used in 
the simulation of Au-t-Au at lA GeV collisions. In one 
previous work, the value 77 = 0.2, i.e. 80% of the free 
space nucleon-nucleon cross section has been obtained in 


B. Effect on radial flow 

Bauer et al. pointed out that in intermediate energy 
HIC, nuclear stopping power is determined by both the 
mean field (EOS) and the in-medium nucleon-nucleon 
cross section [s^. The nuclear stopping was also pro¬ 
posed as a probe to extract information on the isospin 
dependence of the in-medium N-N cross section in HIC 
for the beam energy from the Fermi energy to about 150A 
MeV di. A similar physical quantity, radial flow, is 
studied in this work. Figure shows the px spectra 
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from the Au + Au central collisions at lA GeV with 
different value of the in-medium NNCS factor 77 (0, 0.2, 
0.5, 0.7 and 0.9), fitting well with a function from the 
Blast-wave model. In Figure [SDb), results demonstrate 
that the system has larger radial flow and higher tem¬ 
perature when the in-medium NNCS becomes larger (i.e. 
the lower 77 value). If this scenario is conhrmed and with 
enough precision, then the in-medium NNCS reduction 
factor might be extracted by a comparison with the ex¬ 
perimental result. 




FIG. 7: (Color online) Rep of protons and pions versus pr in 
Au+Au at lA GeV by the IQMD simulation with soft EoS 
and MDI. (a) The Rep of protons w/ and w/o vr absorption 
channel in IQMD model ; (b) The Rep of the charged tt w/ 
and w/o pion absorption channel. 


VI. PION ABSORPTION PROCESS 

As presented in Sec.II, all the pions are produced by 
resonance (A) decay in this IQMD model as the direct 
production is very small in intermediate energy range 
and can be neglected. According to the pion absorp¬ 
tion with a large probability, it will obviously change 
the phas e-sp ace distribution of pions. As discussed in 
Refs.fSTl. l58l|. pion scattering channel influences the an¬ 
gular distribution of pions and the pion absorption chan¬ 


nel plays an important role for the absolute number of 
produced pions. If we turn off this two channels (identify 
with soft-A production channel (ttV —>■ A), here we call 
“pion absorption process”), then pion will emit directly 
from hard-A decay without any scattering and absorp¬ 
tion. The hot medium and cold shadowing effects of pions 
can be quantified by observing Rep. The isotropic angu¬ 
lar distribution of without pion absorption process is 
shown in Figure [ 6 Ka), indicating that the pions are pro¬ 
duced originally via hard A decay, while there exists two 
peaks at 90° and 270° owing to “squeezing out” effect in 
the case with pion absorption. 

In addition, R^p of protons and pions are investigated 
by a comparison between the cases with and without pion 
absorption. In Fig. [3 the upper panel (a) shows the 
Rep of proton w/ and w/o pion absorption in Au-|-Au 
collisions at incident energy lA GeV, respectively; and 
the Rep of changed pion w/ and w/o pion absorption are 
shown in the bottom panel (b). From this figure, we can 
see that, while the Rep of protons has a slight difference 
between the cases w/ and w/o pion absorption, the Rep of 
pion has been changed prodigiously. And the difference 
has demonstrated clearly that “pion absorption process” 
is decisive on pion dynamics. 


VII. SUMMARY 

To summarize, nuclear modification factors Rep of pro¬ 
tons and pions at different incident energies have been 
investigated within the IQMD model with SM-EOS. The 
Rep of protons rises rapidly with the pp increasing at 
0.8A GeV owing to radial flow and Gronin effect. And 
the rising trend becomes slowly with the incident energy 
increasing. This feature can be explained by the nuclear 
stopping whose degree decreases with the incident en¬ 
ergy increase. On the other hand, the Rep of pion rises 
slowly at low pp (< 0.3 GeV) and level off at high pp{> 
O.SGeV). It can be understood by these pions being pro¬ 
duced from two sources, the hard A decayed pions which 
emit outward directly and the soft A decay pions which 
are absorbed and then secondary decayed. The Rep of 
pions has a little enhancement at low pp because these 
low energy pions are affected by nucleon dynamics, such 
as radial flow, and it maintains a saturated trend at high 
Pp owing to no Gronin effect on high energy pions. 

The difference of the nuclear modification factor be¬ 
tween protons and pions indicates the different interac¬ 
tion mechanisms. Because pion has a large absorption 
cross section by nucleon, the pions with high transverse 
momentum are not enhanced in central collisions, but the 
protons are. All these observations are consistent with a 
picture where a dense strongly interacting nuclear mat¬ 
ter with a moderate collective flow is most likely formed 
in central Au -I- Au collisions over a large rapidity range 
which results in the strong suppression of charged pion 
yields at high pp and boosts the protons to higher trans¬ 
verse momentum. 






We change the NNCS in nuclear medium between 0.2- 
0.9. Results demonstrate that radial flow at the central 
collisions decrease with a smaller in-medium NNCS. The 
Rep also became weakly increasing with a lower NNCS. 
Additionally, the Rep obtained in IQMD model without 
pion absorption is investigated in comparison with the 
normal case. The results demonstrate that Rep of pro¬ 
tons have almost no change after deactivating the re¬ 
action channel ttN —>■ A, while the Rep of pions have 


significant difference. This phenomenon reflects that the 
pion absorption mechanism plays a dominat role on pion 
dynamics. 
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